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A random walk lattice model of adsorption of an isolated polymer chain at a solution surface is 
investigated. The model is a modification of a simple cubic lattice in which there is a correlation 
between successive steps. The direction of each step is at right angles to the direction of the preceding 
step (all bond angles are 90°). One-dimensional characteristics of the monomer unit distribution are 
determined analytically in the limit of long polymer chains neglecting the self-excluded volume. The 
mean number of monomer units adsorbed in the surface layer v(0, IV) is determined assuming that one 
end of the polymer chain lies in the surface layer, where N is the mean number of monomer units in 
the chain and is the adsorption energy of each monomer unit in the surface layer measured in units 
of kT. In addition, the mean distance of the free end of the chain from the surface layer z{6, N) is de- 
termined. The properties of this correlated step model are qualitatively similar to the properties which 
have been found in uncorrelated step models. In particular, there is a critical value of the adsorption 
energy C such that for > C , i/(0, N) is proportional to N. Numerical values of N~ l v{6, IV) and z(6, N) 
are presented for > 9 c = \n (V5-1). 

Key Words: Chain polymer, adsorption, random walk, lattice model, short range correlation, 
critical energy, partition function, generating function. 

1. Introduction 

We have recently studied a variety of random walk lattice models of polymer chain adsorption 
at a solution surface [l]. 1 At the same time F. McCrackin [2] has completed a Monte Carlo inves- 
tigation of still another lattice model of polymer chain adsorption at a solution surface in which 
he treats the effect of the self-excluded volume of the polymer chain. The lattice model treated 
by McCrackin is a simple cubic lattice in which the direction of each step must be at right angles 
to the direction of the preceding step. Thus there is a correlation between neighboring steps in 
his model. The purpose of this paper is twofold: {l\ To generalize the methods developed in I in 
order to treat a lattice model with correlation between neighboring steps, and (2) to determine some 
of the moments of the monomer unit distribution which will serve as reference values for the ex- 
cluded volume computations of McCrackin. 

The averages which are calculated are /V ~ x v{6, TV), the average fraction of monomer units ad- 
sorbed in the solution surface, and z(d, N), the mean distance of one end of the chain from the 
surface when the other end is in the surface. The parameter N is the number of monomer units 
in the chain, and is the energy of adsorption of a monomer unit measured in units of kT. The 
lattice model and recurrence equations for the associated generalized random walk are described 
in section 2. The recurrence equations are solved formally in section 3 by the method of generating 
functions; and formulas for N~ l v(6, N) and z($, IV) are obtained in terms of the generating function 
solution. It is shown in section 4 that there is a critical value of the adsorption energy 6 c =\n (V5— 1) 
such that for > 6 C the molecule exists in an adsorbed state. In particular, for 9 > C the fraction 
N~ l v(6, N) is independent of N and greater than zero in the limit TV— » o°. The numerical results 
presented in section 4 for N~ l v{6, N) and z(0, N) cover only the adsorbed state of the chain molecule 
(0 > C ) and apply to the case in which the first two monomer units of one end of the chain lie in 
the solution surface. There are no qualitative differences between the lattice model investigated 
in this paper and the simpler lattice models studied in I. 



i brackets indicate the literature references at the end of this paper. 
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2. Lattice Model and Recurrence Equations 

We consider a simple-cubic lattice model of the solution-surface system in which the solution 
surface corresponds to the x-y lattice plane through the point z = 0. Successive lattice planes 
through z=\, 2, . . . represent the solution phase of the system. Polymer chain configurations 
in the solution correspond to paths generated in a random walk on the lattice between nearest- 
neighbor sites subject to the restriction that the direction of each step must be at right angles to 
the direction of the preceding step. The physical presence of the surface is introduced by con- 
sidering only random walks [3, 4] which never enter the lattice plane through z = — 1. If a polymer 
chain is very far from the solution surface, all random walk configurations of a given length are 
equally likely. However, we are primarily interested in the influence of an adsorbing solution 
surface on the average conformation of a polymer chain. In our lattice model, all random walk 
paths of N steps with n steps lying in the surface layer have the same a priori probability. Relative 
to a random walk configuration of TV steps with n—\ steps lying in the surface, the a priori proba 
bility of a walk with n steps in the surface layer is greater by the factor e e where 6=elkT and e 
is the adsorption energy of a monomer unit. 

For convenience, we use a random walk terminology. Consider the problem of computing, 
for a random walk in which the first two steps lie in the surface layer, the unnormalized or relative 
probability that at the Nth step the random walker is located at lattice point (m x , m y , m z ), where 
m z ^ 0. There are actually three different relative probabilities P x (m x , m y , m z ; 7V+ 1), P y (m x , 
m y , m z ; N+ 1), and P z (m x , m y , m z ; N~\- 1) associated with the three different directions from which 
the random walker arrives at (m x , m y , m z ). These three relative probabilities are related to the 
relative probabilities at the TVth step by the relations 

Px{m x , my, m z , N+ l) = (l/4)(£+ + E^)[P y (m x , m y , m z ; N) + P z {m x , m y , m z ; N)l m z ** 1 (1) 

P y (m x , m y , m z ; N+ l) = (l/4)(£++£-)[P*(m*, m y , m z , N) + P z (m x , m y , m z ; IV)], m z ^l (2) 

Pz{m x , m y , m z , N+ \) = {H%){E+ +E z )[P x {m x , m y , m z ; A0+ P y (m x , m y , m z ; N)l m z ^ 1 (3) 

and 

P x (m x , m y , 0; N+ l) = ^(l/4)(£J + £^)[P 1/ (^, m v , 0; N) + P g (m x , m y , 0; N)] (4) 

P y (m x , m y , 0; 7V+ 1) = e»(l/4)(£+ + E y )[P x {m x , m y , 0; TV) + P z {m x , m y , 0; N)] (5) 

P z (m x , m y , 0; TV+ l) = e (H4,)Et[P x (m x , m y , 0; N) + Py(m X9 m y , 0; N)] (6) 

where E% is an operator which is defined by the relation 

E%Pj(m x , m y , m z ; N)=Pj(m x ±l, m y , m z ; N). (7) 



The operators E% and E^ have bimilar definitions. Equations (1) to (3) describe the relations 
between the relative probabilities at the N+ 1th and the Nth steps outside the surface layer. The 
factor e e in eqs (4) to (6) accounts for the fact that relative to those configurations where m z ^ 1 at 
the 7V+ lth step, the relative probabilities for those configurations where m z = at the /V+lth 
step are greater by the factor e e . The absence of E~ in eq (6) is related to the fact that the random 
walker enters the z = layer from only one direction. 

We are interested in deriving from the solution of eqs (1) to (6) the values of z(6, N) the mean 
distance from the surface at the /Vth step and v{6, N) the mean number of steps in the surface 
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layer. The definitions of these quantities are 



and 



lie, m = & 00, N))lQ(d, N) =j- e In Q(6, N) (8) 



z(6,N)= J J £ J) m z P } {m x ,my, mi ;N)IQ(d,N) (9) 

m x = — x m < .= — x m^=0 j-1 



where 



<?(0,W) = ^ I 2 I Pj<i»x, m„ m,; N) (10) 

m x = — oo m v =-x m 2 =0 j=l 

plays the role of a partition function. The subscript j= 1, 2, 3 denotes respectively, x, y, and z: 
It was noted in 1 and it follows from the structure of the recurrence equations that if one starts 
a random walk with the first two steps in the surface layer and generates in succession the relative 
probabilities P x {m x , m, h m z ; r), P u (m x , m u , m z ; r), and P z (m x , m !h m z ; r), these probabilities are 
linear combinations of exponentials e nd where 2^/i=Sr+2. Differentiation of 0(0, N) with 
respect to has the effect of generating the sum required for computing the average number n 
of steps in the surface layer. 

It should be noted that the probabilities Pj{m x , m lh m z ; N) contain more information than is 
required to calculate the z component of the mean displacement after /V steps and the mean num- 
ber of steps in the surface layer. In fact, we will not solve eqs (1) to (6). Instead we will pro- 
ceed as in I and for each value of m z we will sum the appropriate equations in (1) to (6) over all 
values of m x and m u . The result is 

p x (m z ; N+ l) = (H2)[p u (m z ; lV) + Pz(m z ; N)], m z ^1 (11) 

p fJ (m z ; /V+ l) = (l/2)[p x (m z ; N) + p z (m z ; N)l m z & 1 (12) 

Pz (m z ; N+l)=(U4,)(Ei+Ej){p x (m z ; N)+Py(m z ; /V)], hi* 2*1 (13) 

and 

[l-(l-e-")Jp,.(0; /V+l) = (l/2)[p,,(0; N)+p z (0; N)] (14) 

[l-(l-e-^,(0; N+l) = (U2)\p x (0; N) + p z (Q; N)\ (15) 

[l-(l-^)]p,(0; N+l^a/tyEHpAO; /V) + p,(0; N)] (16) 



^here 



Pj (m z ;N)= £ J P i™x*™v>™z\N). (17) 



w r .= -x m u 



We introduce one further transformation in eqs (11) to (16) 

p(rn; lV) = p z (m z ; N) (18) 

and 

q(m; N) = (U2)[ Px (m g ; N) + p y (m z ; TV)]. (19) 
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In terms of these relative probabilities for being in the lattice plane which is m units from the solu- 
tion surface at the N\\\ step, we can rewrite eqs (11) to (16) as 

q(m; W+l) = (l/2)[q(m; N)+p(m; AD], m 2* 1 (20) 

p(m; N+I) = (ll2)[q(m + l;N) + q(m-I; TV)], m^ 1 (21) 

[l-(l-e-%<0; tf+l)=(l/2)fo(0; N)+p(0; N)] (22) 

[l-(l-e-^)]p(0; 7V-r-l) = (l/2) <7 (l; N). (23) 

Equations (20) and (22) are obtained by combining (11), (12) and (14), (15), respectively. 

The sums in eqs (8) to (10) for v{6, N) and z(#, TV), which involve the Pj(m. r , m lh m z ; 7V)'s, can be 
replaced by simpler sums involving p(m; N) and q(m; N) 

z(6, N)=^ m[2q(m; N) + p(m; N)]/Q(6, N) (24) 

and 

Q(0 9 #) = J [2q(m;N)+p(m; N)]. (25) 

The starting condition for the random walk, namely, that the first two steps are in the surface layer, 
corresponds to the statement 

p(m, 0) = 0, m^O (26) 

and 

q (0 9 0) - e w , q(m 9 0) =0, m ^ 1. (27) 

3. Solution of Recurrence Equations 

The recurrence equations (20) to (23) can be solved by introducing generating functions. 
The procedure used is similar to that used in I. Multiply the equation for q(m, N-\- 1) by (27r)~ 1 f 2 e imf) 
and sum over all values of m ^ 0. The result is 

- (1 - e~°) (27r) 1 / 2 9 (0; 7V+ 1) + 5(0, 7V+ 1) = \ {0> (c/>, N) + <2(0, #)}, (28) 

where 

£(<!>, N) = (2ir) - 1 ' 2 J e^m; W) (29) 

and 

N)=(2tt)- 1 ' 2 ]g e"^p(m; /V). (30) 

»i=o 
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The corresponding equation obtained by combining the equations for the pirn; /V+D's is 

-(1-c-^(2tt)- 1 / 2 M0; /V+l) + ^(<£, N+ 1) =cos <£ 2L($, N) - h {2iT)-^e-^q{^ N). (31) 

Next multiply eqs (30) and (31) by y A + 1 , and sum each over all values of N to obtain the pair 

-(\ -<>-»)(27rr^\k {) ( ))-</((): 0)|+ A (0, v)-^(0, O) = ^y[r(0, y) + A(0, y)] (32) 

and 

-(l-c-^)(277)- 1 / 2 |// (v) -/;((); 0)]+ H0, y)- ^(0, 0) = ycos (f> A((/>, y)-| y^Tr)- 1 ' 2 *?-***:^) (33) 

where 

A(tf>,y)= j? y<2 (<£,#), (34) 

V=0 

n*,y)=J r v ^(<^/V), 05) 

My) =2 >" v /> ( ">; M, (36) 



Ur)= £ fq(m; N). (37) 

.v=o 

The definitions (29), (30), (36), and (37), when combined with the starting conditions (26) and (27), 
are substituted in (32) and (33). The result is 

- \ yn<t>, y) + (l- I y)A(0, y) = (27r)" 1 / 2 [(] -ei%) + ^] (38) 

and 

1(0, y)- vcos 0A(0, y)=(27r)- 1 / 2 [(l -e-^)Ao(y)-|ye-^(y)]. (39) 

Next solve eqs (38) and (39) for T(0, y) and A(0, y) and obtain 

H0, Y) = i x {27T \ '^ r { [( 1 - e- d )Ao( v) - | ye-**A;o(y)]( 1 " \ y) + [( 1 - «>-%,( y) + e*]y cos 0} 

I-2V-2 rcos0 (4Q) 

A(0, v) = - — ; (2?r) 1 ' 1/2 > -r{(l-e- e )kd(y) + e e +%y[a-e- e )ho(y)-h ye^koiy)]}. (41) 

1 — 2 y~2 t cos <p 

Equations (40) and (41) constitute a pair of implicit equations for the functions h {) {y) and k (y) 
because according to the definitions of T((/), y), A(</>, y), h m (y\ and A: m (y) in (34) to (37) and 0» (</>,N) 
and J(0, JV) in (29) and (30) 

F(</>, y) = (27r)- 1 / 2 J) e im *A m (y) (42) 

m = 
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and 



M^y) = (27r)-^ %&*kjy). 



(43) 



A pair of equations for determining A (y) cind A (y) can be obtained by multiplying eqs (40) and 
(41) by (27r)~ 1/2 and integrating with respect to (/> from — it to it 



My) = [d - e- e )h(y)ho(y) - \ yh{y)Hy)]{\ - \ y) + y[(l - e^%) + e%(y) 



and 



(44) 



Arfy) = [( 1 - e-*)&o(y) + e»]/o(y) + 1 y( I - e-*)My)/o(y) - \ y 2 Uy)h{y) 



(45) 



J_ f^ e im(t) d(t) 

Uj) -2vrJ_ 7r . 1 1 2 
1— T^y-^y 2 



y--y 2 cos (/) 



1+I r > + I^ 



-1/2 



9- 7 9 



'-I- 



1+|^+|^ 



1/21 



1 



| m| 



(46) 



Solving (44) and (45) for /i (y) and ko(y), we obtain 



iw[i-< 



i. ,, 



/*<> (y) = j ye»/, (r) | ] — (1 — e-")/ (y) + J y*/,(y) 



ye%iy)h(y) 



|(i-|rl 



/ 0(0, y) (47) 



and 



My)={e"/»(y)[l-(l- |y)(l-e-%,(y)]+L^y)%)(e»-l)}/«(«, y) 



(48) 



where 

9{6, y) = 1 + (1 - e-») (| y- 2W) + J y 2 ?i(j) + d- e-») 2 / (y)[/ (y) (l -| y) - 1 fHy)\ 

Equations (47) to (49) can be simplified with the aid of eq (46) to give 

h n {y)=lyc»I 1 (y)l®(e,y) 
and 

where 



(49) 



(50) 



(51) 



9(6, y) = 



l+\f+\f 



'1 



l+\f+\f 



1/2 11 1 



(52) 
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Substitute the expressions for h,iy) and k (y) obtained in (50) and (51) into eqs (40) and (41) for 
T(4>, y) and A(<f>, y), and then form the quantity V(4>, y)+2A(<j>, y) 



r(<f>, y) + 2A(<f>, y)-- 



(277-)- |/2 (l + iy) 
' [I- i y-$f cos <f> r l* 



|><e 9 -l)/i(y)+/o(y) 



-| y e-^ + 2(l-e-") 2+ 2 y + COS ^ 



9(0, y) 



+ 2e el±£™±\. (53) 

I T V 



Equation (53) constitutes an explicit solution of the recurrence equations (20) to (23) which can be 
used to compute the sums appearing in eqs (24), (25), and (8). In particular, from the definitions 
of T(<fi, y) and A((/>, y) we have 



T(4>, y) + 2A(0, >) = (27t)- 1 / 2 £ r v 2 e h " (l >[p(m; N) + 2q(m; n)]. 



(54) 



A=() Di = 



If we set (/) = () and select the coefficient of y v in eq (54) using Cauchy's formula, we obtain an 
explicit expression for Q(6, N) 



Q(e,N) = -^-.i -^-^TrV^nO, y) + 2A(0, y)] 



(55) 



where C is a closed contour in the complex y-plane which encloses only the pole of the integrand 
at y=0. Substituting eq (53) in (55) and using eqs (46) and (52), one obtains the explicit formula 



1 i dy 1 \hy(eP-\)L(y) + h(y)[2(\-e- 9 )-%y\ 



Q{e ^ = MJ- 



9(6, y) 



+ 2e»\ 



(56) 



Q( , N ) = 2e<>+-Ll ~^~ 

2mJc„y N+i l-y 



(e"-l) 



^+\f+\f 



y-'+2(l-e-")-|y 



1+ I y2 + I r ' 



(57) 



The sum V m[p(m; N) + 2g(m; N)] appearing in eq (24) can be extracted from eq (54) 

| m[p(m; N) + 2q(m; W]=~j c ^l~^j^ ^, y) + 2A(<£, y)] 

The explicit formula obtained from eq (58) after inserting (53) is 

f m[p(m; N) + 2q(m; /V)]=^<f ^ n ^'SL i 
ito 2mJc„y N (l-yW(e,y) 
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(58) 



(59) 



789-972 O - 66 - 4 



or 



£ m[p(m; N) + 2q(m; A)] 



2m J Co y^ 1 — y 



A /(l-y)(l+iy 2 + |y i )-l + ^y-ye-^ + 2e-^ 



(60) 



4. Calculation of *>(0, TV) and z(0, TV) for 7V»1 
In this section we evaluate the contour integrals in eqs (57) and (60) for N > 1 and so obtain 
expressions for v(6, N) from eq (8) and z(0, N) from eq (24). The procedure which we use to evalu- 
ate the contour integrals is similar to that used in I. 

Location of singularities. It is first necessary to locate the singularities of the integrands 
in (57) and (60). In addition to the poles at y=0 and y=l, there are branch point singularities 



associated with zeros of the square root function J(l-y) li-yH-^ 3 ), namely, singularities 
at y— 1 and at the roots of the cubic equation 



4 J 4 J 



(61) 



It is shown in the appendix that these three roots, y„, y/>, and y c , lie outside the unit circle in the 
complex y-plane. In view of the occurrence of these branch point singularities, we introduce 
cuts in the complex y-plane as shown in figure 1 in order to define the integrands in eqs (57) and 
(60) uniquely. The cuts start at the branch points and extend radially outward. The square root 
function is taken as real and positive on the real axis between y = and 1. With this definition 
of the square root function, there may be an additional singularity of the integrands in (57) and (60) 




FIGURE 1. Singularities of integrands, cuts, and contours of inte- 
gration in the complex y-plane. 
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arising from a zero of the denominator 



#(0,y) = A (l-y) 1 + 1^+1^ 



I 



i 



l+-y+2e- 2 "- 



-ye 



(62) 



A numerical investigation of g(6, y) for different values of shows that there is never more than 
one root of g(9< y) inside the unit circle. When there is a root of g(6, y) inside the unit circle, it 
is real, and it lies between y=0 and 1. In order to determine the condition on for which there 
is a real root yi(0) ^ 1 of #(0, y), consider the plots of the two functions F(y) = 



(1 — y) I 1 +Ty 2_ t"Ty 3 ) and L(0, y)= 1 — 2e 2 ^ — (i~ e # )> shown in figure 2. The occurrence 

of the root y\(6) corresponds to an intersection of F(y) and L(0, y). For = ln 2, the line L(0, y) 
becomes L(ln 2, y)= J . For > In 2, the slope of L(0, y) is negative; and as increases the inter- 
section of L(0, y) and F(y) moves up toward the point (0, 1), i.e., toward smaller values of yi(0) than 
yi(ln 2). For < In 2, the slope of L(0, y) is positive; and as decreases the intersection of L(0, y) 
and F(0, y) moves down toward the point (1, 0). The critical value of 0, for which the intersection 
of L(0, y) and F(y) occurs at (1, 0) can be determined by solving the equation g(0, 1) = 0, i.e., 



-l + ! + 2e- 2 "-< 



= 0. 



The value of at which yi(0) = 1 is 



c =ln(V5-l). 



(63) 



(64) 



There are no other singularities of the integrands in (57) and (60). 

Deformation of the contour of integration. Having determined the positions and types of 

singularities of the integrands in the integral representations of (A0, /V) and V m[p(m; N) 

///=o 
+ 2q(m; /V)], we will now evaluate these integrals in the limit in which /V is large and > C . As 
in I, the details of the calculation are different, depending upon whether 0> C , 0< C , or 0=0 C . 



,L(X2,y) 




Figure 2. Plots of F(y)= | (1 - y )(l +\y 2 ^\y 3 )\ and 
L($, y) =\ — 2e~ 29 — (« — e~ e \y versus y for real y, ^ y «S 1. 



309 



The results in each case are qualitatively the same as in I. Therefore, we will limit our discussion 
to the case 6 > 6 C . 

Consider the integral 

1 £ Ay) dy 



J = - 



27TiJg(0,y)y N+1 



(65) 



around the three contours C , C\, and C2 shown in figure 1. Denote these integrals by J r o, J \, and 
J^ 2 , respectively. Suppose that the integrand f(y)lg(6, y)y N+1 is the integrand in eq (57) or (60) 

oc 

so that J() = Q(6, N)or */q = V m[p{m; N)-\-2q(m; M)], respectively. The contour C\ encloses 

yi(#), and the contour C? consists of a circular part whose radius is 1+8, and an indented part 
around the branch point at y= 1. With these definitions, the integrals satisfy the relation 

^0=^2-^/1 (66) 



We first evaluate J \ and indicate why J 2 is negligible compared to J \ for TV >> 1. The value of 
J \ is simply the residue [5] at yi = y\(6) 



Si=fiyi)lg'(0,y l )y«S 



(67) 



where g'(0, y x ) denotes — g(0, y)l y = yi . Thus J \ is proportional to the (N+ l) st power of a number 
ay 

which is larger than one. As in I, it can be verified in the case of the contour integral J 2 that the 

contribution from the circular part of the contour is proportional to (l-h8) _,v , and that the contri- 

bution from the remaining part of the contour is small compared to ( — ) . Therefore, in the 

limit N> 1, we can neglect J 2 compared to J\ and write eq (66) as 



r o = -^i 



(68) 



This approximation is excellent for large finite values of ^V provided that 6 is not too close to d c , 
i.e., provided that y\ is not too close to the branch point y=l. Having established the relations 
(67) and (68) for both integrals (57) and (60), we can now calculate the values of v{0, N) and z(0, N). 
First we have from eq (8) 



fie,N) = ^lnQ(e,N) = -N^/ yi . 



(69) 
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FIGURE 3. Average fraction of adsorbed monomer units 
N _1 v(6, N) versus adsorption energy 6 {solid curve) for > d c . 

For comparison, the corresponding curves for the 6-choice simple cubic lattice (dashed 
curve) and the 2-choice simple cubic (dotted curve) are included, from reference I. 

310 



6.0 




1 


1 1 1 1 1 — 1 


5.0 


- 




- 


4.0 






- 


2 

^"3.0 






\ 


2.0 






\ 


1.0 






v^^ 






< B 


"" — 





1 1 1 1 1 , 



FIGURE 4. Mean distance of end of polymer chain from surface 
z(0, N) versus adsorption energy for 6 > 6 C . 



In eq (69), we have dropped terms which are not proportional to N. It is seen that the result 
obtained here for the case of the 4-choice s.c. lattice when 6 > 6 C is similar to the result obtained 
in I, namely, on the average a finite fraction of all steps are adsorbed in the surface for /V> 1. 
In the case of eq (24) we have 



z(0, N)±- yi \( e °-l) 



1-pv 



i+^+i^ 



yr 1 + 2(l-e-*)-|yi} l , (70) 



a result which is independent of TV f or N > 1 . 

The values of N~ l v(d, N) and z(0, IV) in eqs (69) and (70) have been determined numerically by F. 
McCrackin; and the results are presented graphically in figures 3 and 4 for 6 > 6 C . Values of 
N~ l v{d, N) obtained in I for the 6-choice simple cubic lattice and the 2-choice simple cubic lattice 
are also included in figure 3 for comparison. Although results have only been presented for adsorp- 
tion energies greater than r , examination of the behavior of the integrands in (57) and (60) for 6 < C 
and 6= r in the neighborhood of y=l shows that the results for the 4-choice simple cubic lattice 
are qualitatively similar to the results for the lattices studied in I over the entire range of values of 0. 



I thank F. McCrackin for calculating the numerical values of TV 1 v(6, IV) and z(0, N). 

5. Appendix 

The roots of the cubic equation 



Iy3 + I j2 + 1=r0 
4 J 4 



(Al) 



are the reciprocals of the roots of the equation 



23+72+7=0 

4 4 
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Therefore the roots of (Al) are 
y a = -2{[(l + 1 / 2 7) 1 / 2 +l] 1 / 3 -[(H- 1 /27) 1 / 2 -l] 1 / 3 }- 1 
= -1.968 

y /> = 4{l(l + V27) 1 / 2 + l] 1 / 3 -[(l + V 2 7) 1 / 2 ^ 

= 0.497 + 1.308/ 
y c = 0.497 -1.308* 
Thus all three roots, y (h yt,, and y r , lie outside the unit circle. 
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